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We propose a scheme to deterministically generate Greenberger-Horne-Zeilinger 
states of N > 3 atoms trapped in spatially separated cavities connected by optical 
fibers. The scheme is based on the technique of fractional stimulated Raman adia- 
batic passage which is one-step in the sense that one needs just wait for the desired 



qh. entangled state to be generated in the stationary regime. The parametrized shapes 

of the Rabi frequencies of the classical fields that drive the two end atoms are chosen 
appropriately to realize the scheme. We also show numerically that the proposed 



scheme is insensitive to the fluctuations of the pulses' parameters and, at the same 
time, robust against decoherence caused by the dissipation due to fiber decay. More- 



' over, a relatively high fidelity can be obtained even in the presence of cavity decay 



and atomic spontaneous emission. 
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I. INTRODUCTION 



Entanglement is not on 



y an essential ingredient for testing quantum nonlocality against 



local hidden theories , but also a necessary resource for implementing various quan- 
tum informatic tasks. Fundamentally, it is one of the most important traits in quantum 
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mechanics. It has found different applications in quantum information processing (QIP) 
such as quantum cryptography l3( , quantum teleportation [4j, |5j, quantum dense coding 



7j, quantum secrete sharing [8[, and so on. Typical entangled states are Bell states 



fl, 



Greenberger-Horne-Zeilinger (GHZ) states |2J and W states [9[ have been identified and can 
directly be utilized in QIP [3 8 1 - In particular, great interest has been arisen regarding the 



significant role of GHZ states in the 



and quantum communication 10 
spectroscopy 
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"oundations of quantum mechanics measurement theory 
-|l4|. error correction protocols 15|, |l6j, and high-precision 



181 ] . Contrary to bipartite entangled states, GHZ states exhibit a special 



kind of entanglement between N > 3 parties, providing a possibility to test quantum non- 
locality in a one-shot manner. So, of common interest is the problem of how to gen erate 
GHZ states using current technologies. In fact, for trapped ions [20] and photons [2l|, |22j, 



a series of experimental methods 
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231 ] have already been invented. It is worthy to 



note that cavity quantum electrodynamics (CQED) proves to be very promising for QIP. 
Based on CQED, numerous schemes 24l43l| have been proposed for deterministic genera- 
Ion of entanglement between atoms trapped in different cavities connected by optical fibers 
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36]. For example, Zheng et al. proposed a simplified scheme to product GHZ states 



Both the schemes are 



371 ]. while Li et al. made use of the quantum Zeno dynamics 
difficult to implement because they depend on the exact knowledge of all parameters and 
require controlling the interaction time accurately. A way to overcome such difficulties in 
state engineering is to force the system's initial state evolve along a dark-state, if any, by 
means of adiabatic passage. Such an evolution can be realized by the so-called technique 
of stimulated Raman adiabatic passage (STIRAP) 39| or fractional STIRAP (f-STIRAP) 
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42 1 and 



which have been employed in the context of coherent population transfer 
coherent atomic beam deflection j^ . 

In this paper, we design a one-step scheme to deterministically generate GHZ states of 
N atoms individually trapped in a linear array of optical cavities whose nearest neighbors 
are connected by — 1 optical fibers. Two external lasers are needed to drive the two end 
atom-cavity subsystems. The scheme is based on the adiabatic passage along a dark state, 
which is a specific eigenstate of the total atom-cavity-fiber system corresponding to the zero 
eigenvalue. The key idea is to choose suitable time-dependent Rabi frequencies of the driving 
lasers. Here, we use Gaussian pulses with proper parameters which are turned on and turned 
off in an adequate manner so that in the long-time (stationary) limit the desired GHZ states 
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FIG. 1: Three atoms (a±, 0,2 and 03), each of which has one excited state (\e)k {k = 1>2, 3)) 
and three ground states (\gi)k, \go)k and \git)k) with a tripod-type configuration, are respectively 
trapped in three optical cavities (c%, C2 and C3) connected by two optical fibers (/1 and /2). 



emerge automatically. Compared with the previous schemes [37|, [38| , ours has the following 
advantages: (i) The multi-atom GHZ states are produced deterministically only in one step 
without worrying about precise control over interaction time; (ii) The scheme is insensitive 
to moderate fluctuations of experimental parameters and (iii) The process is immune to the 
fiber decay, and a relatively high fidelity can be obtained even in the presence of cavity 
decay and atomic spontaneous emission. 

Our paper is organized as follows. After the Introduction, in section II, we describe the 
physical model and present the detailed procedure to realize the scheme for generating three- 
atom GHZ states. The general case of any N > 3 atoms is also touched upon briefly in this 
section. Then, in section III, we discuss issues regarding the robustness of our scheme against 
possible fluctuations in the parameters involved as well as against various mechanisms of 
decoherence. Finally, we conclude in section IV. 



II. GENERATION OF MULTI-ATOM GHZ STATES 

We first consider the case of three atoms in detail. As shown in Fig. 1, three atoms 
ax, (I2 and 03 are trapped in three distant linearly arranged optical cavities c%, C2 and C3, 
respectively. Each atom has one excited state \e)k {k = 1, 2, 3) and three ground states \gL)k, 
\go)k and \gR)k, which are nondegenerate and correspond to J = 1 and m = —1, 0, +1, 
respectively. The cavities c\ and C3 are single-mode while the cavity C2 is two-mode. They 
are connected by two short optical fibers fi and f 2 . The atomic transitions |e)i( 2 ) <H- |(?l)i(2) 
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FIG. 2: N atoms are respectively trapped in N cavities through N — 1 fibers. 

(I e )2(3) "H- \gR) 2(3)) are resonantly coupled to the left-circularly (right-circularly) polarized 
cavity mode, while the transitions \e)i -H- \go)i (|e) 3 -h- |^ )s) is driven resonantly by an 
external classical laser. 



In the short-fiber limit, Lu/(2ttc) <C 1 29|, |44( with L the fiber length, c the speed of 



light and u the decay rate of the cavity field into a continuum of fiber modes, only one 
resonant fiber mode interacts with the cavity mode. Then, in the interaction picture, the 
Hamiltonian of the total atom-cavity-fiber system can be written as (h = 1) 

H = H a i + H ac + H c j, (1) 

H al = ^(^leMsol + Q 3 (t)e^ 3 |e> 3 <^o| + H.c, (2) 

2 3 

H ac = '^9i,iaiM) i (g L \ + } j gi, r ai,r\ e )i(9R\ + H.c, (3) 

i=l i=2 

H cf = t>i&l(ai 5 z + a 2 ,i) + v 2 b\(a 2) r + a 3jr ) + H.c, (4) 

with f^i(3)(t) and <fi(3) the Rabi frequencies and phases of the driving lasers, a\ and Oi,i( r ) 
the creation and annihilation operators for the left-circularly (right-circularly) polarized 
mode of cavity q, and (bf) the creation (annihilation) operators of the resonant mode 
of fiber ff. For simplicity, we assume equal atom-cavity and equal cavity- fiber coupling 
strengths, i.e., g it i = g^ r = g and V\ — v 2 — v. We also denote by \n) cl ^ the quantum field 
state of cavity c\ (03) containing n left-circularly (right-circularly) polarized photons, by 
I m, n) c2 the quantum field state of cavity c 2 containing m left-circularly and n right-circularly 
polarized photons, and by /i(/2) t ne quantum field state of fiber fi (f 2 ) containing n 
photons. 

Let the total system be initially in the separable state 

M-OO)) = |^,^,^)a 1 a 2 a3|0)cl|0)/l|0,0) c2 |0) /2 |0) c3 . (5) 
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Then, governed by the Hamiltonian in Eq. (JTJ, it will evolve in a closed subspace spanned 
by 11 basis states {|0z); 1 — 1,2, 11} : 



>i) = I^(-OO)) = |po,^,^)a 1 a 2 a 3 |0) c l|0) / i|0,0) c2 |0) /2 |0) 



c3 • 



102) = \e,g L ,gR) 
10s) 

104) 



|0) c i|0)/i|0,0) c2 |0) /2 |0) C 3, 

\gLi 9l, 9r) aia 2 a 3 |l)cl|0)/l|0, 0) c2 |0)/ 2 |0) c 3, 
\9L) 9l, 9r) |0)cl|l)/l|0, 0) c2 |0)/ 2 |0) c3 , 

105) = |^,^,^)a 1 a 2 a 3 |0) c l|0) /1 |l,0) c2 |0) /2 |0) c3 , 

|0e) = \9L,e,g R ) ai 02(13 

|0) c i|0)/i|0,0) c2 |0) /2 |0) c3 , 
107) = \9l,9r,9r) (11(12(13 

|0) c i|0)/i|0, l) c2 |0) /2 |0) c3 , 

10s) = \9l,9r,9r) £11(12(13 

|0) c i|0)/i|0,0) c2 |l) /2 |0) c3 , 

I0g) = \9l,9r,9r) (11(12(13 

|0) cl |0)/i|0, 0) c2 |0)/ 2 |l) c3 , 
) = \9l, 9r, e) aia2as |0) cl \0) fl \0, 0) c2 |0) /2 |0) c3 , 



and 

I01l) = |^,^^o)a ia2a 3|0)cl|0)/l|0,0) c2 |0) /2 |0) c 3. 

The Hamiltonian ([1]), in terms of the basis states Eqs. ([6] -[16]), reads 
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with n m (t)=n m (t)e-^m 



6 



An alternative basis of the subspace consists of 11 time-dependent eigenstates of the 
instantaneous Hamiltonian given by Eq. ( ITT)) (|$ m (t)) , m = 1,2, 11}. It can be verified 
that one of the eigenvalues of such H is zero and the corresponding eigenstate, which we 
label by |$i(t)) , has the form 

= ° {t) I'M 

y/4X 2 (t) + G 2 (t)[X 2 (t) + iy 

X[t) (103) - 105) + 107) - 109)) 



V4X 2 (t) + G 2 (t)[X 2 (t) + l] 



= l0n), (18) 
'4X 2 (t) + G 2 (i)[X 2 (t) + 1] 

where X(t) = Qi(t)/Q 3 (t) and G(t) = g/Q 3 (t). The state |<3?i(t)) is called a trapped or dark 

state since it contains atoms only in the ground states during the entire evolution. The fact 

that excited levels \e) k are missing in |$i(t)) is due to the destructive quantum interference. 

Namely, as seen from Eqs. (J6] - |T6]) , \4>2) (106), |0io)) does contain the excited level |e) 1 

(|e) 2 , |e) 3 ), but, the transition |0i) — > |0 2 ) (|0s) — > |0e) , 109) |0io)) is canceled by the 

transition \<f) 3 ) — > \<p 2 ) {\4>t) ~~ 106) > |0n) — > |0io))- It is a ls° of a surprise that the fibers' 

modes b\ and b 2 do not appear in |$i(t)) (i.e., |0 4 ) and |0g) that contain a fiber mode are 

absent in This is again a consequence of destructive quantum interference: the 

transitions [fa) — > \<p^) ({fa) —> \4>s)) an d |0s) — >■ (|0g) — > \<f>s)) destroy each other 

completely. However, both the atomic excited levels and the fibers' modes are important, 

via their virtual excitations, in connecting the different atomic ground levels as well as the 

different cavities, a very necessary feature in our system to generate entanglement between 

separated atoms. 

The total atom-cavity- fiber system state at any time t can be extended as a su- 

perposition of {|$ m (t))}, 

n n 

|tt(*)) = X) U>m(t)|$m(*)), ^2\w m (t)\ 2 = l, (19) 
m=l m=l 

with |w m (t)| 2 the probability of finding the system in state \<& m (t)) . If the pulse that drives 
the atom a 3 precedes the pulse that drives the atom ai, i.e., 

lim ft 3 (i) > lim SlAt) = 0, (20) 

t— >— oo t— >— oo 

or lim^_ oc X(t) = X(-oo) = 0, then from Eqs. JTHJ), ([19]) and ([20]) it follows that 
|$i(-oo)) = |0i) = |*(— oo)) . Using this in Eq. (TT9]) yields w m=1 (— oo) = 1 and 
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w ra> i(- oo) = 0. If we vary the pulses slowly enough in time to satisfy the adiabatic following 
condition, then the system initial state should evolve only along of Eq. ffT8]) . The 

key strategy is to tailor the pulses so that 

2n 1 (t)n 3 (t) 

. = V = < g, 21 

(or 2X(t)/ \J X 2 {t) + 1 <C G(t)) to keep the probabilities of finding the system in states 
{|03)> 105)) 107)) 109)} negligible during the evolution and in the long-time limit both the 
pulses vanish simultaneously retaining their ratio constant, i.e., 

lim Q 13 (t) = 0; lim [n 1 (t)/fi 3 (t)] = lim X(t) = X(oo) = const. (22) 



If SO 



with 



lim |$a(t)} = \ghz) 123 ®\0) cl \0) fl \0,0) c2 \0) f2 \0) c3 , (23) 

t— ¥00 



\9hz) 123 = 1 =\gO,9L,gR) ai a 2 a 3 - e 1 ^ 1 ^ ^=^==\g L , 9R, go) ai a 2 a 3 , (24) 

a/a 2 (oo) + 1 a/a 2 (oo) + 1 

which for X(oo) = 1 and ip\ + <y9 3 = rr is the desired GHZ state \GHZ) l23 = (\g , g L , g R ) + 
\9l,9r,9o)) aia 2 a 3 /V2. 

We next briefly present the generalization of the above scheme to the case of iV > 3 
atoms. As Fig. 2 shows, the iV atoms ai,a 2 , <in are respectively trapped in iV cavities 
ci, c 2 , cat connected by iV — 1 fibers fx, f 2 , /jv-i- The level configurations of atoms ax, 
{a 2 , a 3 , cztv-i} and a^ are the same as those of atoms ax, a 2 and a 3 in the case of iV = 3. 
For example, for an odd N > 3 and equal atom-cavity and equal cavity-fiber coupling 
strengths, the expressions of H a i, H ac and H c j in the total Hamiltonian H read 

H al = fi 1 (t)e^|e) 1 ^ | + ^(t)e i<PJV |e)^o| + ff.c., (25) 

PV-1 N 



H ac = g 



22 a iA e )d9L\ + 22 <X iA e )i'\9R\ 



+ H.C, (26) 



. i=l i=2 

{N-l)/2 

H cf = V [&2i-l(°2i-l,J + fl2i,0 + 4i(«2i,r + 02i+l,r)] + H.C.. (27) 

i=l 

Suppose that initially the atoms are prepared in the separable state 
\go, gL, gR, gi, ■ ■ ■ , gR) ai a 2 ...a N while all the cavities and fibers are empty. Then, 
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under the constraint 2Qi (t)il^ (t) / \Z^if (t) + £1% (t) <C g and the adiabatic fol- 
lowing condition, the atoms can eventually appear in the entangled state 
\ghz) aia2 ,„ aN = cosa\g , g L , 9r, 9l, ■ ■ ■ , 9r) + e' l ^ 1+ ^ N hma\g L , g R , g L , g R , . . . , g }) ai a 2 ...a N , if 
lim t ^._ 0O fijv(t) > lim^-oofii^) = 0, Irnif^ f2i,jv(£) = and lim^ 00 [Oi(t)/fi A r(t)] = tana. 



III. REALIZATION AND DISCUSSION 



A possible implementation of our scheme for N = 3 can be realized by using driving 
pulses with the Rabi frequencies of the shapes [40 1 



fii(t) = ^osin« e- (i " T)2/T2 , (28) 

and 

fi 3 (*) = fi e- (t+T)2/T2 + fiocosa e^) 2 /* 2 . (29) 

For illustration, we display in Fig. 3 the time- dependence of the pulses fti^t), the proba- 
bilities P n (t) = | {<fi n | 2 of finding the system in the states \(p n ) (n = 1,3, 11) and 
the fidelity F = \ (GHZ \<&i(t)) | 2 with the parameters a = 7r/4, Q q = O.lg, t = 50/g and 
T = 80/ g. As is clear from the top figure, the pulses (|28[) and (|2"9~1) with the above-chosen 
parameters satisfy the conditions Eq. f )20|) and Eq. (1221) with A(oo) = tana. The middle 
figure shows that Pi is decreasing from 1 to 0.5, Pn is instead increasing from to 0.5, while 
P 3 = P 5 = P 7 = P 9 remains negligible all the time (in fact P3,5 7 g is increasing from until a 
maximum value of about 0.0032 when gt ~ 38 but then is decreasing back to 0). In theory, 
the GHZ state is generated asymptotically in the long-time limit. However, as seen from 
the bottom figure, at gt = 100 the fidelity is already F(gt = 100) ~ 0.99 and at gt = 170 it 
is almost unity: F(gt = 170) ~ 0.9999. 

Next, we will analyze the robustness of adiabaticity condition against the pulse shapes 
of classical fields. From Ref. 40], we know that our scheme needs an optimal range of r 
related to T to achieve a preferable adiabaticity. So during the evolution, a major challenge 
is to choose an optimal relation between r and T. The dependence of the fidelity F on the 
pulses' parameters shown in Fig. 4 indicates that F is larger than 99% for r/T within the 
range 0.37 < r/T < 1.11. This means that our scheme could content with the preferable 
adiabaticity in a relatively large range. 
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0.14 




Time ( in units of g ) 

FIG. 3: The time dependence of (top) the Rabi frequencies of the driving lasers ^l\(t) and £ls(t), 
(middle) the probabilities Pi, P3, P5, P7, P9, P\\ of finding the system in states \<j>\), \<f>s), |</>7), 
l^g), 1 0n) , respectively, and (bottom) the fidelity F. The parameters used are a = 7r/4, f2n/s = 0.1, 
c/r = 50 and c/T = 80. 

So far the whole system is treated as absolutely isolated from the environment, i.e., we 
have totally omitted the decoherence effect in our system. In order to confirm the validity 
of our scheme, we now discuss on the influence of decoherence induced by cavity decay, fiber 
decay and atomic spontaneous emission. To account for the decoherence we resort to the 
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FIG. 4: Density plot of the fidelity F at gt = 300 as a function of gr and gT for flo/g = 0.1 and 
v/g = 10. The two straight lines indicate the boundaries of the region within which the fidelity is 
almost one. 

master equation for the density matrix p(t) of the whole system which has the standard 
form 



P 



[H, P ]-J2^(b f b f p-2b f pb f + pb f b 

2 

E % ( a h a %iP - 2a i,iP a h + P a h a i,i) 

3 



t=2 

-EE f (4^p-2^p4+^)' (3°) 

* =1 j=go,gL,gn 

where iS^. = |e)j(j|, S~j = \j)i(e\, kf (/Cj) denotes the decay rate of the fibers (cavities) 
and 7y is the spontaneous emission rate of the atoms. We assume kf = k, Ki = k and 
7y = 7 = 7o/3 for simplicity. The fidelity F versus the ratios flo/g and n/g (fto/ 9 an d 1 /g) 
is displayed in the top (bottom) panel of Fig. 5. We can see from Fig. 5 that with the 
increasing of the laser intensity, the decoherence caused by the atomic spontaneous emission 
is getting smaller and smaller, while the decoherence caused by the cavity decay is becoming 
greater and greater. The reason is that the adiabatic passage is just likely to evolve within 
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FIG. 5: Density plot of the fidelity F at gt = 300 as a function of (top) &o/g and n/g and (bottom) 
tt /g and 7/5. 

the dark state subspace under relatively large laser intensity, and to violate the influence 
caused by the spontaneous emission as well. However, the probabilities that the cavity fields 
are excited increase with laser intensity, which in turn increase the dissipation caused by the 
cavity decay finally. An appropriate value Qq should be chosen by taking into account both 
the factors (spontaneous emission of atoms and decay of cavities) as the two error sources 
cannot be avoided simultaneously. The fidelity F versus the cavity decay n/g and the fiber 
decay k/g is shown in Fig. 6, where we neglect the spontaneous emission of atoms. As 
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FIG. 6: The fidelity F as a function of n/g and k/g for Vt^/g = 0.1, u/g = 10 and gt = 300. 




0.05 



FIG. 7: The fidelity F as a function of n/g and 7/5 for Qo/g = 0.1, u/g = 10 and gt = 300. 

seen from the figure, the fidelity F decreases with the increasing cavity decay, but is almost 
unaffected by the fiber decay. Even though we set k/g = 0.1 (and 7/5 = n/g = 0) the 
fidelity is still as high as F = 0.998, so the decoherence due to the fibers hardly influences 
the quality of the generated state. The fidelity F versus the decay of cavities n/g and the 
spontaneous emission of atoms 7/gr is shown in Fig. 7 with the fiber decay ignored. We see 
from Fig. 7 that F decreases with the increasing of both the cavity decay and the atomic 
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spontaneous emission. For a relative large value of 7/(7 = n/g = 0.05, the fidelity is still 
about F = 0.850 when the laser intensity Qo/g — 0.1 chosen in our scheme. Therefore, our 
scheme is robust in realistic conditions. 

Finally, let us discuss on the experimental feasibility. The parameters g = 2tt x 75 
MHz, 7 = 2tt x 2.62 MHz and k = 2tt x 3.5 MHz are achievable in optical cavities with 



the wavelength in the region 630 — 850 nm in recent experiments [45J, |46| . A near-perfect 
fiber-cavity coupling with an efficiency larger than 97.20% can be realized using fiber-taper 



coupling to high-Q si 
is about 2.2 dB/km 



ica microspheres 



471 ]. The optical fiber decay at a 852 nm wavelength 



48l . |49j], which corresponds to the fiber decay rate of k = 0.152 MHz, 
which is lower than the cavity decay rate. With these parameters, we will obtain a high 
fidelity F, meaning that it is possible to realize our scheme in a realistic experiment. 



IV. CONCLUSION 



In conclusion, we have proposed a one-step scheme for deterministic generation of GHZ 
states for any number of atoms individually trapped in spatially separated cavities connected 
with short optical fibers via adiabatic passage by appropriately tailoring the external driving 
laser fields. The figure of merit is that we need not to precisely control the generation time: 
the desired state emerges as a steady state of the system. It is interesting that the generation 
time (i.e., the time it takes to generate the GHZ state with a desired quality, i.e., F = 1 — e 
with a predetermined small e) does not increase with the number of atoms. This is of 
importance from the view point of decoherence when dealing with a large-sized system 
such as entanglement of a big number of atoms. We have also numerically calculated the 
influences of the driving laser's parameters as well as of the decoherence effect caused by the 
atom spontaneous emission and the cavity/fiber decay on the quality of the generated state 
in terms of fidelity. The numerical results have revealed that a relatively high fidelity of 
three-atom GHZ states can be obtained in the presence of the above-mentioned influencing 
factors. Therefore, we do hope that within the current experimental technology it would be 
possible to realize our scheme. 
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